4.2 Fisher的信任推断法

Fisher 的区间估计的方法原则上可以用于任何统计推断问题, 代表了对统计问题的一种根本上不同的观点.

1 信任分布

设样本 X∼N(θ,1), 样本大小为 1, 则 X−θ∼N(0,1), 即对 ∀t∈R: P(X−θ<t)=Φ(t)=12π∫−∞te−y22dy, 它可以改写为 (1.1)P(θ>X−t)=Φ(t),或P(θ<X−t)=1−Φ(t). 虽然从通常的概率论来说这样的写法没什么区别, 但是 Fisher 却认为, 有了样本 X 后, 把 θ 看成一个随机变量, 它就有了一个分布, 把它称为 θ 的信任分布.

在得到 X 前, 我们对 θ 一无所知. X 提供的信息打破了我们的无知, 用概率分布的形式给出了对 θ 的新认识.

再比如, 要做 θ 的 区间估计, 给定 α, 找 a,b(a<b), 使 P~(a≤θ≤b)=1−α (P~ 表示 θ 的信任分布). 则根据 (1.1), 这要求 Φ(X−a)−Φ(X−b)=1−α. 取一组让 b−a 最小的 a,b, 这时 a=X−uα2, b=X+uα2.


对于多个样本 X1,⋯,Xn∼i.i.dN(θ,1), 注意到 X― 是 θ 的 充分统计量, 且 n(X―−θ)∼N(0,1), 则得到信任分布 θ∼N(X―,1n). 由此建立信任区间 [X―−uα2n,X―+uα2n].
再比如 X1,⋯,Xn∼i.i.dN(a,σ2), a,σ 未知, 则 n(X―−a)S∼tn−1, 即 P(n(X―−a)S<x)=Tn−1(x)=P(a>X―−Sxn).
虽然到目前为止和置信区间的结果一样, 但是往后人们发现具体结果也可能不同.
Fisher 的理论目前存在两个问题:

2 用 Fisher 方法解 Behrens-Fisher 问题

Behrens-Fisher 问题是这样的问题:

问题

X1,⋯,Xm∼N(a,σ12), Y1,⋯,Yn∼N(b,σ22), a,b,σ1,σ2 未知, 所有样本独立, 找 b−a 的区间估计.

也即和前面的两样本检验相比, 最大的不同是方差不再相同.
为了解决这个问题, 记 X―,S12,Y―,S22. 记 t1∼tm−1, t2∼tn−1. 用 ξ=dη 表示两个随机变量分布相同. 则 m(X―−a)S1=dt1,n(Y―−b)S2=dt2. 记 Z=Y―−X―, θ=b−a. 记 S1∗=S1m, S2∗=S2n, 则 Z−θ=S2∗t2−S1∗t1⇒θ=Z−(S2∗t2−S1∗t1). 有了样本后, 记 Z,S1∗,S2∗ 的具体值为 z,s1∗,s2∗. 则上式变为 θ 的信任分布.
因为 t1,t2 独立且都是 t 分布, 所以可以确定 θ 的信任分布. 记 r=s1∗2+s2∗2, 找 ψ:cos⁡ψ=s2∗r. 此时 sin⁡ψ=s1∗r. 而 s2∗t2−s1∗t1=r(t2cos⁡ψ−t1sin⁡ψ). 用 Fm,n,ψ 表示 t2cos⁡ψ−t1sin⁡ψ 的分布(因为只和 m,n,ψ 有关). 找 ym,n,ψ,α>0 使得 Fm,n,ψ(ym,n,ψ,α)−Fm,n,ψ(−ym,n,ψ,α)=1−α, 则 P~(|z−θ|≤rym,n,ψ,α)=P(|t2cos⁡ψ−t1sin⁡ψ|≤ym,n,ψ,α)=1−α, 于是 P~(z−rym,n,ψ,α≤θ≤z+rym,n,ψ,α)=1−α. 从而我们得到信任区间 [z−rym,n,ψ,α,z+rym,n,ψ,α].
可以证明这个区间并不是 1−α 置信区间.